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Abstract. The dark matter content of galaxies is usually determined from galaxies in dynamical 
equilibrium, mainly from rotationally supported galactic components. Such determinations restrict 
measurements to special regions in galaxies, e.g. the galactic plane(s), whereas other regions are 
not probed at all. Interacting galaxies offer an alternative, because extended tidal tails often probe 
outer or off-plane regions of galaxies. However, these systems are neither in dynamical equilibrium 
nor simple, because they are composed of two or more galaxies, by this increasing the associated 
parameter space. We present our genetic algorithm based modeling tool which allows to investigate 
the extended parameter space of interacting galaxies. From these studies, we derive the dynamical 
history of (well observed) galaxies. Among other parameters we constrain the dark matter content 
of the involved galaxies. We demonstrate the applicability of this strategy with examples ranging 
from stellar streams around the Milky Way to extended tidal tails, from proto-typical binary galaxies 
(like M5 1 or the Antennae system) to small group of galaxies. 
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INTRODUCTION 

Interacting galaxies often show extended arms and bridges, caused by the gravitational 
tides [1], [2], [3]. Modeling these interactions allows us to recover the kinematic history 
of observed galaxies and serves as a test bed for the physics of the merging process. 
Timescales, dynamical friction, galactic structures and many more can be investigated 
and verified by models. By observing rotation curves of spiral galaxies, Zwicky [4] 
has shown, that galaxies must contain an additional non visible mass. This dark matter 
(DM) was introduced to galaxy simulations mostly as a collision-less particle compo- 
nent. Shape and extension of the (DM-)halos have large effects on the behavior of in- 
teracting galaxies [5]. State-of-the-art simulations include multiple components (stars, 
DM, different gas components, magnetic fields, radiation) and complex recipes for star 
formation and feedback, AGN influence etc. These simulations require a huge amount 
of computational power. For remodeling observed systems it is therefore important to 
know as many physical parameters of the interaction as possible in advance. To solve 
this problem, we coupled a fast restricted N-body code for galactic interactions to a 
genetic algorithm (GA). 
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FIGURE 1. Left plot: Interaction of equal disk galaxies. Deflection from the parabolic orbit (for point 
masses) depending on (isothermal) halo size. From right to left rhaio = Okpc to rhaio = 140 kpc in steps of 
20kpc (centers marked with diamonds). Right plot: same as left, but including dynamical friction. 



METHOD 

Finding interaction parameters means solving a minimization problem in an high dimen- 
sional parameter space. We use a GA based on PIKAIA [6] and an improved restricted 
N-body code based on Toomre & Toomre [3]. One important step of improvement was 
the implementation of static (DM-)halos (e.g. isothermal, NFW or Burkert profiles) in- 
stead of point masses for galaxies. The evolution of our galaxies in time was changed 
from solving the Keplerian orbit to a time dependent evolution of (overlapping) poten- 
tials. Test particles mimicking the disk are introduced and serve as imprints for the tidal 
features compared with observed galaxy's HI features. Recently, we have implemented 
dynamical friction to the code to allow for later stages of interaction and merging [7]. 

STRUCTURAL PARAMETERS OF DM HALOES 



Our code was already successfully applied to derive tidal models for NGC4449 [8], 
M51 [9] an the Magellanic Clouds [10]. Here we want to focus on the latter and also on 
the effects of the DM halo and dynamical friction. 

The size (and shape) of the DM halo has an important influence to the deflection from 
a Keplerian orbit and tidal features of interacting disk galaxies. By comparing observed 
features with our models, we can derive these halo parameters. Fig. 1 (left) shows the 
large change in deflection angle and tidal features depending on the halo extension, 
where the final distance keeps almost unaffected. The latter changes when dynamical 
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friction (accounting for density gradients, mass and distance) is switched on (Fig. 1, 
right). 



The Magellanic Stream and the DM halo of the Milky Way 

The most recent proper motion measurements for the Large (LMC) and the Small 
(SMC) Magellanic Clouds [11, 12] put them on significantly larger galactocentric ve- 
locities. We confirmed the failure of the tidal stripping models of the Milky Way 
(MW)-LMC-SMC interaction in such a case [10]. That happened due to the interaction 
timescale reduced as the Clouds' orbital periods become comparable to the Hubble time. 
Shattow et al. [13] have resolved the orbital issue by increasing the MW circular velocity. 
It affected both the LMC/SMC spatial motion and the MW total mass, leading to mul- 
tiple perigalactic passages for the Clouds. We have analyzed the parameter space of the 
Magellanic interaction by a robust search algorithm (OA) combined with a fast 3 D N- 
body model of the tidal interaction involving a flattened dark matter halo of the Galaxy. 
The LSR circular velocity was varied between 210km s~^ and 260 kms^^ and linked 
with the LSR galactocentric distance by the LSR angular rotation rate of 29.45 km s^^ 
[14]. Under such assumptions the Clouds were able to stay confined in the MW halo 
virial radius for well over 4 Gyr, and the neutral hydrogen (HI) Magellanic Stream [15] 
was reproduced quite well over the entire range of the LSR circular velocities. The tidal 




FIGURE 2. Left plot: Orbital evolution of the Clouds over the last 4 Gyr for a high-quality model with 

Vsun = 211.2kms"\ q = 0.74, Mmw = 0.7 ■ lO'^Mj . Right plot: Visualization of the entire simulated 
three-dimensional H I data cube of the Magellanic System. The column density isosurface E = 5 • 
ax is shown, together with the data cube projected to the 2 D map of the integrated column density 
in the position-radial velocity space. The dashed line in the position-velocity projections depicts the LSR 
radial velocity gradient along the observed Magellanic Stream. 

model worked with the new LMC/SMC proper motions over the MW halo mass range 
of O.5-31O^^M0. In agreement with Ruzicka et al. [16], an oblate (flattening<1.0) log- 
arithmic halo was preferred (see Fig. 3). However, the modeled Magellanic Stream was 
constantly displaced in its projected position compared to the observations. Introducing 
the LSR angular rotation rate according to [14] increases the Galaxy mass and reduces 
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the LMC/SMC galactocentric velocities in their magnitude, and so the tidal model comes 
back to play. But it also alters the LMC/SMC velocity vectors in their direction causing 
an undesired shift of the Stream in the position-LSR radial velocity space. 

The studied parameter space of the Magellanic interaction involves the velocity of 
the Clouds (/iN, Hw, Vrad)^ LMC/SMC particle disk radii, galactic circular velocity 
at the solar position (Vsun)> and the MW halo flattening (q). With the use of GA, 

10^ parameter combinations were tested in total, and 100 sets providing satisfactory 
reproduction of the H I Magellanic large-scale structures were collected. 
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FIGURE 3. Distribution of 100 GA fits of the Magellanic System over the analyzed ranges of the 
flattening q of the logarithmic dark matter halo of the Galaxy, and of the solar galactocentric circular 
velocity Vq™. The plot also shows the isolines of the total mass of the Milky Way dark matter halo 
enclosed within the radius of 250 kpc. The mass of the logarithmic dark matter halo is a function of 3 
independent parameters out of which the flattening q and the solar galactocentric circular velocity Vq™ 
were considered free parameters. 
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